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Comment on “Chain Length Scaling of 
Protein Folding Time” 

In a recent Letter, Gutin, Abkevich, and Shakhnovich 

[1] reported on a series of dynamical Monte Carlo sim- 
ulations on lattice models of proteins. Based on these 
highly simplified models, they found that four different 
potential energies lead to four different folding time scales 
Tf, where r/ scales with chain length as N x (see, also, 
Refs. [2-4]), with Å ranging from 2.7 to 6.0. However, 
because of the lack of microscopic models of protein fold- 
ing dynamics, the interpretation and origin of the data 
have remained somewhat speculative. It is the purpose of 
this Comment to point out that the application of a sim- 
ple “mesoscopic” model of protein folding [2] provides a 
quantitative interpretation of the data, as well as a major 
qualitative difference with Ref. [1]. 

The main features of the theory [2] are as follows. Con- 
sistent with the heuristic arguments regarding the en- 
tropic cost of loop formation [1,3], we model the acqui- 
sition of native-like features by considering the full set 
of loops or links of an ideal chain of N monomers, and 
a target State defined by a unique set of N/2 “native- 
like” links. Folding is tuned by a microscopic frustration 
parameter A, defined as the contact energy ratio 8/s be- 
tween random and native-like links, respectively. For a 
finite chain, the model has an effective folding transition 
temperature Tf(N) ~ (1 — A)/lnJV, between random 
conformations and native-like structures. At Tf(N), a 
bimodal distribution on the nurnber of non-native con- 
tacts is observed, whereas away from Tf(N) distributions 
are unimodal. We point out that the latter is commonly 
associated with “cooperativity” [1] . However, we predict 
the size of the critical folding nucleus to be proportional 
to the size of the system. The same conclusion is found 
based on exact enumerations of self-avoiding walks [5]. 
In agreement with the ideal “minimally frustrated” Go 
model X Go ~ 2.7 [1], optimal folding to fully equilibrated 
native-like structures takes place in an entropic domi- 
nated time scale Tf N x ° , where Åo — 3. The striking 
prediction of the theory is that at the transition Tf(N) 
— which, in all likelihood, corresponds to the “optimal” 
temperature estimated using rnean first passage time sim- 
ulations [1] — T f (Tf) ~ N x , where Å ~ 4 (see Fig. 1), 
the same behavior now being confirmed by simulations 
[1] on specially designed fast folding sequences! 

This scaling breaks down when the frustration limited 
folding time scale ta ~ N x ° exp(asA/T) becomes larger 
than Tf, a ~ 0(1) being a constant. These scaling forms 
predict that, for small enough A < 1/3, proteins fold fast 
in a unique range of temperatures roughly independent 
of size and structural specificities, i.e. Ta (IV) < T < 
Tf(N), where Ta is a dynamically determined transition 
temperature at which Tf ~ ta- Below Ta the relaxation 
is slow, following an Arrhenius law [1] (higher energy 
barriers not included here may also be important in this 
regime). If Tf(N) < Ta (IV) frustration dominates alto- 
gether, and folding can be extremely slow [2], with 


Tf (Tf) ~ i\r Ao+2A / (1-A ). (1) 

Simulations have shown that different potential ener- 
gies have different landscapes with different degrees of 
frustration [1,4]. The model averages frustration into a 
single model dependent parameter A. As sketched in Fig. 
3 of Ref. [2] , random or frustration dominated sequences 
with A $ 1/2 will not be able to fold in a physically rele- 
vant time scale. Equating \ RAN ~ 6 with the exponent 
in (1), we find /\ RAN ~ 3/5, which is quite appropriate 
for a random selection of contact energies (cf. [2]). 

It has been claimed that these results can be under- 
stood in the context of a first-order-like folding transition 
[1] . However, it is well known that near a first-order tran- 
sition different relaxation times grow in different ways. 
In particular, these authors have also argued about the 
existence of a macroscopic barrier between native and 
non-native States (“cooperativity”). Thus, if this bar- 
rier needs to be crossed over, the relaxation must then 
be exponential in IV! On the contrary, we find an effec- 
tive second order folding transition, where at Tf(N) any 
typical relaxation time must grow as a power of N. In 
summary, the data in [1] is accounted for by [2], giving 
further credence to the cross-linking model, and empha- 
sizing the role of loop formation in folding kinetics. 
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FIG. 1. Folding time scale at the transition as a function 
of size N. The slopes correspond to 3.8 ± .25 (A = 0) and 
4.3 ± .2 (A = 1/3) [2]. For details of master equation see [2]. 
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